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We examine the quality of the local self-energy approximation, applied here to models of multiple 
quantum impurities coupled to an electronic bath. The local self-energy is obtained by solving a 
single-impurity Anderson model in an effective medium that is determined self-consistently, simi¬ 
lar to the dynamical mean-field theory (DMFT) for correlated lattice systems. By comparing to 
exact results obtained using the numerical renormalization group, we determine situations where 
“impurity-DMFT” is able to capture the physics of highly inhomogeneous systems, and those cases 
where it fails. For two magnetic impurities separated in real-space, the onset of the dilute limit is 
captured, but RKKY-dominated inter-impurity singlet formation cannot be described. For parallel 
quantum dot devices, impurity-DMFT succeeds in capturing underscreened Kondo physics by self- 
consistent generation of a critical pseudogapped effective medium. However, the quantum phase 
transition between high- and low-spin states on tuning interdot coupling cannot be described. 


I. INTRODUCTION 

Dynamical mean field theory (DMFT) has been well- 
established as a powerful tool for treating correlated 
materials.^ At its heart is a local self-energy approxima¬ 
tion that only becomes exact in the limit of infinite lattice 
coordination.^ In many applications to real materials,^ 
direct comparisons with experiment have shown the qual¬ 
ity of this approximation to be very good. However, the 
interacting lattice fermion models to which DMFT is ap¬ 
plied typically have no known exact solution with which 
to fully compare. 

Although DMFT was originally formulated for 
translationally-invariant problems, the local self-energy 
approximation can also be used in an inhomogeneous 
context. Each inequivalent site in the system is mapped 
to a single-impurity Anderson model in an effective 
medium that is determined self-consistently.^ This in¬ 
homogeneous (or real-space) DMFT has been applied 
to a diverse range of systems, including materials with 
surfaces^, layered heterostructures,atoms in optical 
traps,and even nanostructures or devices for molecu¬ 
lar electronics. 

Inhomogeneous DMFT has also been used to study 
depleted lattices.In the case of metals doped with 
randomly distributed magnetic impurities, the physics is 
dominated by effects of disorder; while heavy fermion 
physics is recovered on approaching the periodic limit, 
where standard DMFT applies. In the other limit of di¬ 
lute impurities, inhomogeneous DMFT is exact by con¬ 
struction since the impurity effective medium is simply 
the clean host metal. 

The physics of dilute magnetic impurities is of course 
also well-known,with the impurity spin degrees of free¬ 
dom being screened by surrounding conduction electrons 
by the Kondo effect at low temperatures.^^ Models of 
single magnetic impurities can be solved exactly with the 
numerical renormalization group^^^^^ (NRG), which pro¬ 
vides access to thermodynamic and dynamic quantities 
on essentially any temperature or energy scale. 


However, the physics of the intermediate regime with 
a few coupled quantum impurities is notoriously com¬ 
plex. Already in the two-impurity case, there is a subtle 
competition between single-impurity Kondo physics and 
through-host RKKY coupling.Theoretically, multi¬ 
impurity problems are inherently much more difficult 
to treat: A'-impurity models generally involve N cou¬ 
pled conduction electron channels. For NRG, the re¬ 
quired computational resources scale exponentially with 
the number of channels. In fact, full dynamical proper¬ 
ties of the two-impurity system have only recently been 
obtained with NRG for a true real-space system, tak¬ 
ing into account details of the host lattice exactly.Al¬ 
though the efficiency of NRG can be greatly enhanced 
by interleaving different Wilson chains^^’^^ or exploiting 
large symmetries,^^ it is clear that an exact solution for 
true many-impurity systems is beyond reach. Similar dif¬ 
ficulties afflict other numerical methods. 

Goupled quantum dot devices^^ and magnetic nanos¬ 
tructures or single-molecule junctions^^”^^ can also be 
described theoretically in terms of generalized quantum 
impurity problems. However, these systems typically 
contain many correlated ‘impurity’ degrees of freedom 
(and several effective screening channels), again essen¬ 
tially precluding an exact numerical treatment. 

Inhomogeneous DMFT is therefore a highly appealing 
approximate technique for such systems, since the com¬ 
putational complexity grows only linearly in the number 
of inequivalent interacting sites.^ However, details of ge¬ 
ometry (or molecular structure), inter-impurity coupling 
and RKKY effects are important in these contexts. It is 
therefore not at all obvious what the quality of the local 
self-energy approximation is, since these systems are far 
from an infinite-coordination limit. 

In this paper we apply inhomogeneous DMFT to two 
paradigmatic quantum impurity problems — two mag¬ 
netic impurities separated in real-space, and systems of 
coupled quantum dots between metallic leads (see Fig. 1 
and Sec. II). As with regular DMFT, the local self¬ 
energy scheme in this “impurity- DMFT” represents a 
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rather sophisticated approximation (Sec. Ill), and its ap¬ 
plication here produces nontrivial results. Even though 
the real-space two-impurity model and parallel quantum 
dot models are among the simplest examples for which 
DMFT can be used, they are in fact already at the limit 
of complexity for an exact treatment with NRG. These 
models therefore present an almost unique opportunity 
to test the validity of the local self-energy approximation 
by comparing exact NRG and impurity-DMFT results. 

In the context of quantum impurity models. Ref. 29 
focused on the two-impurity Anderson model in a one¬ 
dimensional geometry. Calculation of various static 
spin susceptibilities (and comparison with DMRG re¬ 
sults) showed that the impurity-DMFT is quantitatively 
correct at large impurity separations. However, the 
impurity-DMFT clearly breaks down at shorter dis¬ 
tances, when it predicts a transition to a symmetry- 
broken state. 

Here we focus on dynamical quantities; in particular 
spectral functions, which are the central experimental ob¬ 
servables for these systems. In the two-impurity model, 
we consider explicitly a three-dimensional (3D) cubic lat¬ 
tice host. The exact NRG solution discussed in Sec. IV A 
shows a range of correlated electron physics, with RKKY- 
dominated effects for small impurity separations, and 
standard Kondo in the dilute impurity limit. Impurity- 
DMFT is correctly able to determine the onset of the di¬ 
lute limit,and certain low-energy properties; however 
when the Kondo effect is suppressed by non-local inter¬ 
impurity RKKY coupling, impurity-DMFT fails quali¬ 
tatively (see Sec. IV B). On the other hand, impurity- 
DMFT is shown to perform surprisingly well for models 
of parallel coupled quantum dots, capturing for example 
nontrivial underscreened Kondo physics through the self- 
consistent generation of a critical pseudogapped effective 
medium in the local problem — see Sec. V. Magneti¬ 
zation effects on application of a dot field are also very 
well captured (Sec. VC), although the quantum phase 
transition between high- and low-spin dot states on tun¬ 
ing inter-dot coupling is not describable within impurity- 
DMFT (Sec. VD). We conclude in Sec. VI with a wider 
discussion of the applicability of inhomogeneous DMFT. 


II. QUANTUM IMPURITY PROBLEMS AND 
EXACT SOLUTION WITH NRG 


Quantum impurity systems, involving either one or 
several coupled impurities, are modeled in terms of gen¬ 
eralized Kondo or Anderson Hamiltonians.^^ They can 
be decomposed as = i^host + ^imp + ^hyb, where 
i^host and i^imp describe the isolated host and impurities, 
respectively; and i^hyb describes the coupling between 
them. In the following, we focus on models of N coupled 
Anderson impurities. Each impurity a is regarded as a 
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FIG. 1. Schematic of the models and mappings used, (a) 
Two separated impurities on the surface of a 3D cubic lattice, 
mapped to a problem involving two Wilson chains in NRG. (b) 
Mapping to a single-impurity problem in a self-consistently 
determined effective medium, via the local self-energy approx¬ 
imation (“impurity-DMFT”). (c) Quantum dots coupled in 
parallel between source and drain leads, mapped to a multi¬ 
impurity but single-channel problem in NRG. 


single spinful correlated quantum level. 
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where creates an electron on level a = 1,..., V with 
spin a The single-particle level energy is 

while Ua is the local Goulomb repulsion penalizing double 
occupancy. Additionally, there may be inter-impurity 
couplings, embodied by idimp-imp- In the case of coupled 
quantum dots. 


d^imp-imp — '^a/3^acr^/3cr T ’ (^) 

a^l3,a 


where v^is is a tunnel-coupling matrix element connect¬ 
ing different dots a and /3. For real magnetic impurity 
clusters on surfaces, Eq. 2 describes a direct intra-cluster 
coupling. In modeling multilevel quantum dots, the dif¬ 
ferent levels are regarded as different ‘impurities’; then 
^imp-imp may include a ferromagnetic Hund’s rule ex¬ 
change between levels.Gapacitive couplings between 
different dots could also be included. 

The host metal is taken to be noninteracting; in a di¬ 
agonal basis it is described generically by, 

^host = ^ ^k^\cT^ko- • (^) 

k,a 
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In the thermodynamic limit, the conduction electrons 
form a continuous band of half-width and total den¬ 
sity of states per spin ptot(<^) = ~ 

The hybridization between the impurities and conduc¬ 
tion electrons is given by, 

= + (4) 

OC k,(T 

Note that the corresponding Kondo models are ob¬ 
tained via a Schrieffer-Wolff transformation^^ in the limit 
of single occupation for each level, perturbatively elim¬ 
inating excitations involving empty or doubly-occupied 
levels to second order in i^hyb- 

The question of how many electronic screening chan¬ 
nels are involved in such models is a subtle one. Although 
in general one can expect Nc = N screening channels for 
a system of N impurities, in fact any 1 < Nc < N can 
be realized, depending on the parameters W/c- 

Consider the N x N impurity Green function matrix 
which contains all the information on the im¬ 
purity single-particle dynamics. In general, there are fi¬ 
nite local and non-local elements = ((d^cr? 

(as usual, {{A;B))^ is the Fourier transform of the re¬ 
tarded correlator The matrix 

Dyson equation then reads, 

[G™P(c^)] = [g™P(^)] . (5) 

Here, X)cr(<^) is the interaction self-energy matrix, con¬ 
taining all of the nontrivial correlated electron physics 
of the problem. The noninteracting (but host-coupled) 
impurity Green functions are given by 

[gyP(a;)] ■' = + iO+)I - E - V - r{oj) , (6) 

where [E\ap = Ma/3 = '^a/3 5 and elements 

of the hybridization matrix are given by = 

In certain high-symmetry cases, the hybridization ma¬ 
trix r(co’) can be diagonalized by a single canonical trans¬ 
formation of operators (that is, independently of energy 
uj). A given channel in this basis is then decoupled 
from the impurity subsystem if it has zero eigenvalue, 
tii{uj) = 0, yielding Nc < N. 

For Nc = N one generally expects an overall singlet 
ground state, with each impurity eventually fully Kondo 
screened. By contrast, for Nc < N there is the possibil¬ 
ity of underscreening, in which the ground state remains 
degenerate. We consider examples of each below; see also 
Fig. 1. 


We will assume here that there is no direct inter¬ 
impurity coupling, idimp-imp = 0. We take as a concrete 
physical example, a semi-infinite 3D cubic tight-binding 
host with a ( 100 ) surface, 

= E(4<.S- + H-c.), (7) 

(b‘>,cr 

where t is the tunneling matrix element connecting 
nearest-neighbor sites i and j. This simplified real-space 
realization captures many features of real metals — in 
particular, there is a definite bandwidth 2D = 12t, and 
the surface density of states is finite and fiat around the 
Fermi level, with p\oc{^ = 0) = l/( 6 t). Eq. 3 is the 
diagonal representation of Eq. 7. 

In this real-space basis, the hybridization Hamiltonian, 
Eq. 4, reads, 

77hyb ^ = E E + H.C. . (8) 

Ck: = l ,2 cr 


Elements of the hybridization matrix then follow as 
ra/?(w) = V*Vi^G^{v„,yi 3 ,uj), where G°(r„,r/ 3 ,w) = 
is the propagator between sites Fq, and 
of the clean host (without impurities). Due to trans¬ 
lational invariance on the surface of the clean host, 
rii(w) = r22(w) = rioc{w) and ri2(w) = r2i(w). 
Note that the non-local functions depend on 

the inter-impurity separation R = ri — r 2 , and must 
be computed explicitly for each real-space realization of 
the two-impurity problem. Here we compute the lattice 
Green functions numerically exactly using the convolu¬ 
tion method of Ref. 19, which can be formulated for sys¬ 
tems with a surface. 

We take now equivalent impurities, = e, Ua = U, 
and Va = V. By transforming to an even/odd impurity 
orbital basis dcjo,a = :^(<^icr i <^ 2 ^)? fhe 2 x 2 matrix 
Dyson equation, Eq. 5, is diagonalized to yield. 


gimp 

'^e/o,o 



= w-hf0+ -e-re/o(w) -Se/o.a(‘^) 


(9) 


where ftc/o{^) = ^ii(<^) A ^ 12 (<^) for ft = G, T, and E. 
The even and odd channels are strictly decoupled in the 
non-interacting {U = 0) system: the even(odd) impurity 
combination couples only to the even(odd) conduction 
electron channel. However, note that Ee(co’) and Eo(co’) 
each contains information on both impurities, and their 
mutual correlations. 

Since re(co’) 7 ^ ro(co’) 7 ^ 0 , the full real-space two- 
impurity Anderson model ( 2 IAM) is irreducibly two- 
channel; see Eig. 1(a). 


A. Two-impurity problem in real-space 

When two Anderson impurities a = 1 and 2 are em¬ 
bedded in a host metal at different real-space positions 
ri and r 2 , a rich range of correlated electron physics can 
result due to the competition between the Kondo effect 
and a through-host RKKY coupling. 


B. Coupled quantum dots 

Several coupled quantum dots, connected between 
source and drain metallic leads, are modeled by Eqs. 1-4. 
Here we focus on the specific geometry of quantum dots 
in paralleP^ [see Eig. 1(c)], which also mimics a multilevel 
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dot.^^ We consider explicitly the symmetric setup with 
Ca = Ua = U^ and Vij = v. The dot-lead hybridization 
for a system of N dots is taken to be, 

N 

<yb'' = ^EE4<.Co. + H.c., ( 10 ) 

a.=l (7 

where Cq^ is an operator for the symmetric combination 
of lead orbitals coupling to the dots. The other conduc¬ 
tion electron states couple indirectly to the dots via Cq^. 

The impurity Green functions for this parallel quantum 
dot (PQD) model are given by the N x N matrix Dyson 
equation, Eq. 5, with elements of the hybridization ma¬ 
trix in Eq. 6 identical by symmetry, = rioc(<^)- 

Eor simplicity, in the following we take rioc(<^) to be the 
same as in Sec. II A. 

In fact, only the symmetric (‘even’) combination of dot 
orbitals couples directly to the leads, 

= VNxVJ 2 dlco. + H.c. , (11) 

a 

with deer = ^cicr- The other dot orbital 

combinations couple indirectly through dea- The ma¬ 
trix Dyson equation can again be diagonalized, yield¬ 
ing in particular for the even-orbital Green function 

G™P = ((4<,;dU)c., 

[G)(y(a;)]“' =w + z0+-e-(7V-l)?;-re(w)-Se,<.(w) , 

( 12 ) 

where re(co’) = A^rioc(<^)- The other {N — 1) con¬ 
duction electron channels decouple from the impurities, 
ro(co’) = raa{(jo) — rc^^/ 3 (co’) = 0, leaving a single-channel 
description. Note, however, that the self-energy Ee(co’) = 
^aa{^) + still inherently contains inter¬ 

impurity correlations. 

C. Numerical renormalization group 

A wide range of quantum impurity problems can be 
solved exactly using NRG,^^ including the above real- 
space 21 AM and PQD models. 

The NRG method involves mapping the full quantum 
impurity model to a ID form amenable to iterative diag- 
onalization. Specifically, the impurities must be coupled 
to the end of one or more decoupled Wilson chains — see 
Eig. 1. This necessitates a diagonal structure of the hy¬ 
bridization matrix Tagioo) oc da/s, each nonzero element 
of which is discretized logarithmically, and mapped to 
a semi-infinite ID tight-binding Wilson chain.Start¬ 
ing from the impurity subsystem, Wilson orbitals are 
successively coupled on, and the Hamiltonian defined in 
the larger Hilbert space is diagonalized. Only the Nx 
lowest energy states are retained for construction of the 
Hamiltonian at the next step. This truncation scheme is 
justified by the exponential decay of hoppings down the 
Wilson chain. Successively lower energies are reached at 
each step, revealing the RG character of the problem. 


In the case of the real-space 2IAM, r(co’) is diagonal 
in the even/odd orbital basis, with = rii(co’) ± 

ri 2 (co’), such that information about the real-space sep¬ 
aration of the impurities is encoded in the difference be¬ 
tween re(co’) and ro(cc;). Discretizing each and mapping 
to even and odd Wilson chains leads to the setup illus¬ 
trated in Eig. 1(a), with two impurities and two channels. 
In our NRG implementation, we discretize using 

A = 2.5, label states by total charge and spin projection 
quantum numbers to block diagonalize the Hamiltonian, 
and retain Nx = 15000 states at each step. So-called 
^-averaging was not required. 

Generalized two-channel problems of this type are 
computationally demanding.Eurthermore, the local 
environments of the two impurities must be identical to 
accomplish the above mapping. If, for example, the im¬ 
purities are different distances away from the surface, or 
if the host is disordered, then r(co’) cannot simply be 
diagonalized by canonical transformations. Similarly, in 
the case of three impurities, NRG can only treat a sym¬ 
metric geometry with = rioc(<^) and ra^g{(jo) = 

rnon-ioc(<^)- The resulting three-channel problem would 
also require significantly greater computational resources 
than that of the two-impurity case. 

In the case of PQD models, the formulation in Sec. H B 
is single-channel: re(co’) can be discretized and mapped 
to a Wilson chain with the dots at one end, as illustrated 
in Eig. 1(c). The complexity here is contained in the 
impurity Hamiltonian itself, which involves N interact¬ 
ing orbitals. In practice this limits NRG to models with 
N < S. The possibility of Kondo underscreening and 
frustration due to finite interdot coupling v means that 
a rather large number of states Nx = 5000 — 8000 must 
be kept at each step to accurately capture the low-energy 
physics (for A = 2 ). 

Impurity spectral functions presented in the following 
sections are obtained from the matrix Dyson equation, 
Eq. 5, using the appropriate exact non-interacting propa¬ 
gators gJf^P(co’). The nontrivial physics enters through the 
N X N self-energy matrix, which is calculated here using 
NRG via a generalization of Ref. 33, from the full den¬ 
sity matrix^^ constructed in the complete Anders-Schiller 
basis. 


III. LOCAL SELF-ENERGY APPROXIMATION 
FOR QUANTUM IMPURITY PROBLEMS 

While the above models can still be solved exactly 
by standard NRG, they are computationally demand¬ 
ing, and represent an approximate limit on what can 
be achieved at present. Glearly true many-impurity sys¬ 
tems and complex nanostructures are out of reach us¬ 
ing full NRG. This provides motivation for development 
of approximate schemes capable of dealing with Kondo 
physics in multi-impurity systems. 

An inhomogeneous (or real-space) DMET approach 
has been widely used to treat generalized Hubbard mod¬ 
els which lack translational invariance.^ At its heart, it 
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involves a local self-energy approximation, which allows 
the full problem to be mapped onto a set of auxiliary 
single-impurity Anderson models (one for each inequiva¬ 
lent site in the system), each in an effective medium that 
must be determined self-consistently. Quantum impurity 
problems can in fact be regarded as highly inhomoge¬ 
neous Hubbard models, and we apply a similar local self¬ 
energy approximation (or “impurity-DMFT”) here to the 
real-space 21 AM and PQD models — see Fig. 1 . 

We first recap briefly the generalized method. The 

impurity Green functions are obtained from the matrix 
Dyson equation Eq. 5 using the exact non-interacting 
propagators gJ^P(co’). However, within impurity-DMFT 
the self-energy matrix takes the form, 

. (13) 

For each inequivalent site 7 , we introduce an auxiliary 

Anderson impurity model, 

"" ~ ( 14 ) 

+ E SvLcfca + E(bfc4cfe^ + H.C.) , 

k,a k,a 

where the parameters and are specified by the 

original quantum impurity model, e^k and V^k are how¬ 
ever determined self-consistently. All the relevant infor¬ 
mation about the effective medium to which the impu¬ 
rity couples is contained in the hybridization function 
Aj(cj) = |t^/cp/(co’ + iO+— e^/c). The impurity Green 

function of this auxiliary model is therefore given by, 

[7T(^)]“' =uj + iO+-e^- A,(c.) - • (15) 

The self-energy is then calculated (e.g. using 

standard single-site single-channel NRG) for each in¬ 
equivalent site 7 , given the hybridization function Aj(cj). 

Aj(cj) is determined self-consistently on the level of 
the local Green functions by demanding: 

, (16a) 

= I] 7 (^) • (16b) 

The self-consistency loop is as follows: 

( 1 ) Guess an initial form for Aj(cj). 

( 2 ) Galculate for auxiliary model Eq. 14, using 

Aj(cj) as input. Repeat for all inequivalent sites 7 . 

(3) Equate E!:^^(c(;) = construct E)cr(co’) from 

Eq. 13, and then find Gl^^^(cu) from Eq. 5. 

(4) Equate and find a new A^{uj) 

via Eq. 15. 

(5) Repeat steps (2)-(4) until converged (i.e. GI^p^(co’) 
does not change upon further iteration). 

An important feature of this self-consistency, which 
is perhaps not immediately obvious, is that it depends 
crucially on the basis used for Eq. 5. 

Eor example, in the real-space 2 IAM, different solu¬ 
tions can be converged, depending on whether one uses 


G-™^(cc;) with i and j labeling either the physical impu¬ 
rity sites 1 and 2 , or their even and odd orbital com¬ 
binations. Indeed, in the latter case, there is no self- 
consistency loop at all: comparing Eqs. 9 and 15 yields 
directly Ae/o(co’) = foi* fhe inequivalent even and 

odd impurities. The local Green functions for the phys¬ 
ical impurities (which are identical by symmetry) follow 
simply as the average, G\™^(cc;) = 4 [G^J^p(co’)+ GJ,^p(co’)]. 

By contrast, establishing self-consistency directly on 
the level of G\™^(co’) implies that Ai(co’) = rii(co’) + 
[ri 2 (cj)]^/[co’ + — e — rii(cL;) — E\^^(cc;)]. Iterating the 

DMET cycle then produces nontrivial structure in Ai(co’) 
due to feedback from the self-energy. This is the formu¬ 
lation employed in Sec. IV. 

Similarly, for the PQD models studied using impurity- 
DMET in Sec. V, we establish self-consistency for the 
local Green functions in the physical impurity basis. 

In this work, we use NRG^^ to solve the auxiliary An¬ 
derson model, Eq. 14. Eor all impurity-DMET calcula¬ 
tions presented, we use A = 2 and Nk = 3000. 

IV. REAL-SPACE TWO-IMPURITY MODEL 

We examine now the real-space 21 AM, with two impu¬ 
rities located at different sites on the (100) surface of a 3D 
cubic lattice [Eqs. 1, 7 and 8 with N = 2; see Eig. 1(a)]. 
Eor simplicity we focus on the particle-hole symmetric 
case, e = —1//2. We recap first the true physics of 
the model, with exact numerical results obtained from 
the full NRG solution, before comparing with impurity- 
DMET results. Eig. 2 shows a comparison for the impu¬ 
rity spectral functions, with exact results as solid black 
lines, and impurity-DMET as red dashed lines. 

A. Full NRG solution 

At very large inter-impurity separations |R| ^ (X), 
each of the two impurities behaves independently and 
one recovers the well-known Kondo physics of the single¬ 
impurity Anderson model.Indeed, the ‘dilute limit’ of 
essentially independent Kondo screening sets in rather 
rapidly with increasing |R| in 3D systems, as known 
from experiments^ and theory.In the dilute limit, the 
impurities are each screened by surrounding conduction 
electrons below an effective single-impurity Kondo tem¬ 
perature At low energies and temperatures, the 

impurity spectral resonance (experimentally observable 
in scanning tunneling spectroscopys^) takes a universal 
Kondo form.s^ 

For two impurities on the (100) surface of the 3D cubic 
lattice, the dilute limit is effectively realized for an inter¬ 
impurity separation R = {Rx^Ry) = (0,4) (with lattice 
constant ao = 1). The bottom-right panel of Fig. 2 shows 
the corresponding impurity spectral function, which is of 
single-impurity form. Note in particular that the Friedel 
sum rule^^ (FSR) implies spectral pinning in the sin¬ 
gle impurity case —Im riocG^^™P(co’ = 0 ) = 1 , which as 
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FIG. 2. Comparison between exact NRG results (solid black 
lines) and impurity-DMFT results (dashed red lines) for the 
real-space 21 AM. The local impurity spectral function at T = 
0 is plotted as —Im FiocG\’^^(ct;) u/D. The two impurities 
are taken to be on the (100) surface of a (semi-) infinite 3D cu¬ 
bic lattice, separated by the vector R = {Rx,Ry). Clockwise 
from top-left panel: R = (0,1), (1,1), (0,4), (2, 2). Impurity 
parameters: U/D = 0.5, e = —U/2 and V/D = 0.075. 


such provides a measure of the independence of Kondo 
screening in problems with several impurities. Here, 
Floe = —Im rioc(O). For R = (0,4), the FSR is very well 
satisfied by NRG results for the full two-impurity prob¬ 
lem. Furthermore, the non-local propagator is 

negligible for this separation (not shown), indicating van¬ 
ishing inter-impurity correlations. 

On bringing the impurities closer, however, genuine 
two-impurity effects manifest due to through-host RKKY 
correlations.^^ This has a marked effect on the impurity 
spectral functions, as seen clearly in the other panels of 
Fig. 2. In addition to differences in the detailed form 
of the impurity spectrum, note that its Fermi-level value 
depends on R, and is not pinned by the FSR. Indeed, the 
nontrivial behavior of — Im riocG\™^(0) is well-known^^ 
from simplified models of the 21 AM; its origin is two- 
impurity contributions to the self-energy 5]^(cc;). 

For separations R = (1,1) and (2,2), the effective 
RKKY coupling is ferromagnetic; the impurities bind to¬ 
gether into a triplet state on the scale of | Jrkky |, which 
is ultimately exactly Kondo-screened by the even and 
odd conduction electron channels. The RKKY scale can 
be read off from Fig. 2 as the energy scale where the 
spectrum departs from single-impurity behavior. As ex¬ 
pected, Jrkky is smaller in magnitude for R = (2, 2) 
than (1,1). The anisotropy between even and odd chan¬ 
nels in the case of R = (1,1) is rather significant (due to 
large Fi 2 (cj)), leading to effective two-stage quenching of 
the inter-impurity triplet. For R = (2,2), the even/odd 
anisotropy is weaker, and two-channel Kondo screening 
occurs essentially in a single step. 



uj/D for the systems plotted in Fig. 2. 


For neighboring impurities R = (0,1), the effective 
RKKY interaction is antiferromagnetic. The impurities 
bind together into a singlet state on the scale of | Jrkky |, 
and decouple from the host, giving G\^^(0) ^ 0. 

Note that in all cases, the overall ground state is a 
spin-singlet (there is no phase transition on varying pa¬ 
rameters or impurity separation^^). The lowest-energy 
physics is Fermi-liquid like, with ^ uj‘^ as 

cj ^ 0. 


B. Impurity-DMFT results 

We turn now to our impurity-DMFT results for the 
same systems — see red dashed lines in Fig. 2. 

In the dilute limit, one expects effective single-impurity 
physics. The full self-energy matrix 5]cr(<^) then becomes 
diagonal, meaning purely local. The impurity-DMFT 
scheme is exact in this limit. This is in fact seen in 
the bottom-right panel of Fig. 2, for impurity separa¬ 
tion R = (0,4): the exact NRG and impurity-DMFT 
solutions are perfectly in agreement. This is confirmed 
in Fig. 3, where the DMFT effective medium A(co’) is 
plotted. For R = (0,4), we find under DMFT self- 
consistency that A(co’) Fioc(co’), meaning that the pres¬ 
ence of the second impurity does not affect the local elec¬ 
tronic environment of the first. Although this might seem 
trivial, it should be noted that the impurity-DMFT is 
correctly able to identify the onset of the dilute limit 
at finite (and surprisingly small) impurity separations. 
Recovering single-impurity physics in impurity-DMFT is 
only strictly trivial as |R| ^ oo. 

The impurity-DMFT scheme is however seen to break 
down as the impurities are brought closer together, and 
inter-impurity effects (embodied by non-local self-energy 
contributions) become important. The failure becomes 
increasingly severe on decreasing the separation. 

For R = (2,2) and (1,1), the spectral behavior for 
energies |c(;| ^ | Jrkky | is correctly described within 
impurity-DMFT; but the effective two-channel screening 
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of the inter-impurity triplet is not captured at lower en¬ 
ergies. In particular, the two-stage screening in the case 
of R = (1,1) is not apparent in the impurity-DMFT re¬ 
sults, with the effect that the Kondo temperature (char¬ 
acterizing RG flow to the Fermi liquid fixed point) is sig¬ 
nificantly overestimated. Tk is modified from the pure 
single-impurity result, but not sufficiently to correctly 
simulate the exact two-impurity results. 

However, it should be noted that in all Kondo-screened 
cases, the impurity-DMFT correctly recovers the nontriv¬ 
ial Fermi level value of the spectrum. This is achieved 
under DMFT self-consistency by the generation of a non¬ 
trivial effective medium, with A(co’ = 0) deviating from 
rioc(<^ = 0)7 a.s demonstrated in Fig. 3. This is con¬ 
sistent with Ref. 29, where an impurity-DMFT scheme 
was used to study two impurities in a ID chain. In that 
work, static quantities at T = 0 were considered, and 
shown to reproduce DMRG results for large separations. 
In 3D, as considered here, the dilute limit sets in much 
more rapidly. These results suggest that in the Kondo 
regime, quantities that depend only on G\^^(0), such as 
the T = 0 resistivity due to impurities, might still be 
accessible with impurity-DMFT. 

Finally, for neighboring impurities R = (0,1), we see 
a catastrophic failure of impurity-DMFT. Here, non¬ 
local self-energy contributions are responsible for inter¬ 
impurity singlet formation, and cannot be captured by 
impurity-DMFT. In fact, impurity-DMFT predicts a 
Kondo-like resonance, with enhanced (rather than sup¬ 
pressed) spectral intensity at the Fermi energy. 


V. COUPLED QUANTUM DOTS 



FIG. 4. Exact NRG results for the PQD model consisting of 
N coupled quantum dots. Plotted is the (even orbital) im¬ 
purity spectral function —Im riocGe’^(ct;) vs uj/D at T = 0 
for N = 1,2, 3,4, 5 (black, red, blue, green and purple lines, 
following the arrow). Inset: rescaled spectra satisfying pin¬ 
ning condition, Eq. 19. Impurity parameters: U/D — 0.5, 
e = —U/2, V/D — 0.06 (and u = h = 0). 


which (in the case of particle-hole symmetry) implies. 


N = 1 


Im -{w/Tk? , 


Im F 


loc 




( / rji \2 

- {uj/Tk) , 


Im riocG™P,(0) = 1 . 


(17a) 

(17b) 

(17c) 


We now discuss systems of N quantum dots, coupled in 
parallel between two metallic leads, Eqs. 1-3 and 10 — see 
Fig. 1(c). Exact NRG results are surveyed in Fig. 4, while 
a detailed comparison between exact NRG and impurity- 
DMFT results is presented in Fig. 5. In the following 
we use a comparatively large U/D = 0.5 to allow for 
a detailed analysis of the universal low-energy physics. 
Although this parameter regime would be appropriate to 
describe the Kondo physics of e.g. 3d transition-metal 
impurities in break junctions, real quantum dot devices 
are typically characterized by a smaller U/D. We note 
that impurity-DMFT can also be applied to the latter 
case with equivalent accuracy. 


A. Full NRG solution 

For a single quantum dot between source and drain 
leads {N = 1), the standard Kondo effect is known to 
arise,producing a dot spectral function that takes 
a universal form when rescaled in terms of T/Tk and 
uj/Tk- The dot spin-^ degree of freedom in this case is 
exactly screened by the lead conduction electrons below 
Tx, forming a ground state spin-singlet. The low-energy 
physics is described in terms of Fermi liquid theory, 


The FSR is embodied by Eq. 17c. This, and the charac¬ 
teristic quadratic approach of the spectrum to its Fermi 
level value Eq. 17b, can be seen from exact NRG results 
in Fig. 4 as the black line. 

However, the physics of PQDs with N > 1 (and simi¬ 
larly, multilevel dots^^’^^) is rather different.In par¬ 
ticular, an 0{V^) effective RKKY interaction arises be¬ 
tween the dots, mediated by the conduction electrons, 

-^RKKY ^ -rlju. (18) 

Importantly, this coupling is ferromagnetic: on temper¬ 
ature/energy scales ^ I JrkkyI, the N dots form a com¬ 
bined spin-A/2 object. This spin-A/2 is coupled to a 
single channel of conduction electrons, and so is under¬ 
screened to a spin-(A — l)/2 object by the Kondo effect 
below Tk ^ Trkky- The ground state therefore remains 
spin-degenerate^^ for A > 1. The conduction electrons 
still feel a 7r/2 phase shift (at T = 0 and u; = 0) due 
to the Kondo effect, resulting in a generalized pinning 
condition for the even orbital spectral function, 

-Im AriocG^7(0) = 1 . (19) 

This is demonstrated in Fig. 4; the inset shows that 
Eq. 19 is satisfied numerically at low energies for all A. 
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FIG. 5. Comparison of impurity-DMFT results (red lines) and exact NRG (black lines) for PQD models with N = 2,3,4, 5. 
Main panel in each case shows the rescaled even orbital spectrum —Im FiocGe™c 7 ^(ct;) vs cj/Tk at T = 0, while the lower insets 
show the unsealed (raw) spectra. Upper insets show the detailed low-energy behavior. Here Tk is defined via Im = 

Tk) = 0.95 X Im = 0). Dot parameters: IJjD — 0.5, e = —Ul2, V/D = 0.05 (and u = h = 0). 


However, the full scaling form of the dot spectrum does 
depend on — see Fig. 4. In particular, there are sin¬ 
gular corrections^^ at low energies due to residual fer¬ 
romagnetic correlations between the underscreened dot 
spin-(A^ — 1)/2 and conduction electrons for \uj\ C Tk- 
This results in a ‘marginal Fermi liquid’ at low energies, 

N >1: 

Im Ee,<.(w) “ - hr^{uj/TK) , (20a) 

Im r,oc - UT(0)] “ ■ (20b) 

The low-energy spectral behavior is thus rather different 
for A/" = 1 compared with N > reflecting the difference 
between exactly screened and underscreened Kondo. 

B. Impurity-DMFT results 

In the case of a single dot N = impurity-DMFT 
is exact by construction. In this subsection, we explore 
the extent to which the local self-energy approximation 
in impurity-DMFT can capture nontrivial underscreened 
Kondo physics in PQDs for N > 1. Fig. 5 compares ex¬ 
act NRG results (black lines) with impurity-DMFT re¬ 


sults (red) for the even orbital spectral function in PQD 
models with A^ = 2, 3,4, 5. 

The impurity-DMFT solution is certainly seen to be 
approximate: the value of Tk is always underestimated 
(lower insets of Fig. 5), and the precise scaling form is 
not entirely recovered (main panels). However, impurity- 
DMFT is apparently rather successful at capturing cer¬ 
tain aspects of PQD physics. 

At higher energies |ci;| ^ Tk, one observes 

1/ln^(co’/Tx) spectral tails as expected. Furthermore, 
for |u;| ^ Jrkky, this behavior should be characteristic 
of single-impurity scattering. Perhaps it is not surpris¬ 
ing the impurity-DMFT captures this physics asymptot¬ 
ically, since the PQD is mapped to an effective single¬ 
impurity model. 

Far more interesting is the low-energy physics, which 
impurity-DMFT emulates surprisingly well in all cases 
considered (see upper insets of Fig. 5 for each N). This 
is a rather nontrivial fact, since the exact results display 
marginal Fermi liquid physics^^ whose physical origin^^ 
is the Kondo underscreening of a high-spin object. This 
type of Kondo effect is obviously not accessible within a 
single impurity Anderson model. 

Instead, the low-energy physics of the PQD model is 
mimicked within impurity-DMFT through generation of 
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FIG. 6. Top panel: DMFT effective medium vs 

uj/D for PQDs with N — 2, 3,...,10 increasing in the direction 
of the arrow. Low-energy behavior follows a pseudogap pow- 
erlaw A(ct;) ~ with exponent a ~ N~^ that depends on 
the number of dots (see inset, points). We find numerically 
7 ~ 0.43 (inset, dashed line). The bare hybridization is shown 
as the red dotted line for comparison. Bottom panel: corre¬ 
sponding local DMFT self-energy —Im Ej^^(ct;), showing the 
same low-energy power law decay as A(a;). Same parameters 
as Fig. 5. 


an effective medium with nontrivial low-energy structure. 
This is shown in the upper panel of Fig. 6 for different 
numbers of coupled dots, N = 2,3, ...,10 (increasing in 
the direction of the arrow). Although at high energies 
|<^| ^ the DMFT effective medium A(ci;) ^ Fioc(<^) 
is essentially that of the bare hybridization (dotted line), 
at lower energies Ic^l < Tk a pseudogap develops in A(u;). 
Specifically, we find the low-energy behavior to be well- 
described by, 

AH |:^r, (21) 

where the exponent a = a{N) > 0 depends on the num¬ 
ber of dots. Although there may be logarithmic correc¬ 
tions to this pure powerlaw behavior, we could not es¬ 
tablish that from the numerical results. As shown in the 
inset to the upper panel of Fig. 6 , the exponent itself is 
found to follow a powerlaw form, a ^ N~^ with 7 ^ 0.43. 
In particular, we note that 0 < a < ^ for all A" > 1. 

It is well-known that the particle-hole symmetric pseu¬ 
dogap Anderson model with 0 < a < ^ supports a 
quantum phase transition separating Kondo strong cou¬ 
pling and local moment phases.The critical point 
itself is interacting and is characterized by nontrivial 
correlations; the signature of this is that the impu¬ 
rity self-energy vanishes with exactly the same power 
as the hybridization.Interestingly, we find within 


impurity-DMFT for PQD models that, 

-Im . (22) 

with the same exponent a as in Eq. 21. This is shown in 
the lower panel of Fig. 6 , and confirms that under self- 
consistency, impurity-DMFT generates a critical pseu- 
dogapped effective single-impurity model. The resulting 
even orbital spectral function"^^ then appears to give a 
very good approximation to the low-energy behavior aris¬ 
ing in underscreened PQD models — see Fig. 5. 

The appearance of an effective pseudogapped local 
model for PQDs within impurity-DMFT is rather remi¬ 
niscent of the results of Ref. 45. In that work, a parallel 
double dot system was considered, with one interacting 
and one non-interacting dot. An exact local model was 
derived, integrating out the non-interacting conduction 
electron bath and also the other non-interacting dot; the 
effective hybridization was found to be pseudogapped 
with exponent a = 2. Our results therefore generalize 
these findings to fully interacting PQDs, albeit that our 
local models arise within the (approximate) impurity- 
DMFT scheme. Since impurity-DMFT is exact by con¬ 
struction for the case considered in Ref. 45, it would be 
interesting explore this connection by tracking the evolu¬ 
tion of the pseudogap as interactions on the second dot 
are switched on. 


C. Magnetic field 

We consider now the effect of applying a magnetic 
field h to the dots. The ‘impurity’ Hamiltonian Eq. 1 
is supplemented by the term iLmag = where 

= |Ea[4tG “ is the total PQD spin 

projection. This Zeeman term can be straightforwardly 
incorporated into the existing formalism by simply re¬ 
placing e ^ ecr in Eqs. 12 and 15, where Ca- = c — haj2 
(and cr = ± for t / i)- 

In the standard single-impurity case N — 1, an in¬ 
finitesimal field h does not destroy the Kondo singlet 
(which has a binding energy ^ Tk)^ and the impurity 
remains unmagnetized. On the other hand, for h ^ Tk 
spin-flip scattering, which is at the heart of the Kondo 
effect, is suppressed. 

In the present context of PQDs with A > 1 , the physics 
is more subtle: even for infinitesimal /i, the residual dot 
spin-(A — 1 )/ 2 , resulting from Kondo under screening, be¬ 
comes polarized. In terms of the dot spectral function, 
this implies a significant (spin-dependent) redistribution 
of spectral weight between the Hubbard satellites — an 
effect not observed in the single-impurity case. The same 
basic physics has been discussed for multi-level dots in 
Ref. 46. Eor larger h ^ T^, the Kondo effect is again 
destroyed. Most interesting is the case where the applied 
field h ^ Tx, so that the Kondo effect competes with spin 
polarization. Non-trivial physics associated with partial 
Kondo screening is then expected. 
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N = 2 



FIG. 7. Comparison of impurity-DMFT results (red dashed 
lines) and exact NRG (solid black lines) for the PQD model 
with N = 2 (upper panels) and 3 (lower panels) in an applied 
magnetic held. Dot parameters U/D = 0.5, e = —f//2 and 
V/D = 0.06 such that ~ 10~^D. Shown for u = 0 but 

h/D = 10“^°, 10“^, 10“® and 10“^. Plotted for a =t- 


In Fig. 7 we present exact NRG results (solid black 
lines) for PQDs with N = 2 (upper panels) and N = 
3 (lower panels), in the case of hnite held hjTK ^ 
0.1,1,10,100. This is the intermediate regime where the 
Kondo effect is only partially manifest due to competing 
spin polarization. In all cases, spectral weight is shifted 
between Hubbard satellites, indicative of dot magnetiza¬ 
tion. As h increases, the Kondo resonance is progres¬ 
sively diminished, and a held-induced peak develops at 
uj ^ h. The dot held cuts ohF the logarithmic corrections, 
and regularizes the spectrum for |cc;| <^h. 

The impurity-DMFT results are shown in Fig. 7 as the 
red dashed lines, and compare extremely well. Several 
nontrivial features of the full N-dot model are recovered 
by impurity-DMFT. 

In particular, the redistribution of spectral weight to 



FIG. 8. DMFT ehective medium ys uj/D for 

PQDs with A = 2, 3,4 and 5 (increasing with the arrow) for 
U/D = 0.5, e = -H/2, V/D = 0.06 and h/D = 10“® {v = 0). 


the upper Hubbard satellite is captured essentially ex¬ 
actly. This might seem surprising, because impurity- 
DMFT involves the mapping to a single-impurity model 
in which small helds are not expected to produce a 
large magnetization. Instead, the effect emerges from 
self-consistent development of structures in the effective 
medium A(u;) on the scale of |ci;| ^ 7/, as shown in Fig. 8. 
The true PQD magnetization is in fact captured very well 
in impurity-DMFT. 

The Fermi level value of the spectrum in Fig. 7 is also 
recovered with high accuracy by impurity-DMFT. This is 
a nontrivial result of the renormalized effective medium 
at the Fermi level, A(cc; = 0) — see Fig. 8. 

In addition, the spectral peak dXuj ^ h (arising in these 
underscreened PQD models due to the asymmetric spin- 
resolved Kondo resonance"^^) is described very well by 
impurity-DMFT. Fig. 8 again shows that this is a result 
of the self-consistent feedback of the local self-energy into 
the effective medium. 

Indeed, the entire energy dependence of the spectrum 
is very well described by impurity-DMFT for both N = 2 
and N = 3 over the full range of fields considered. The 
quality of the results improve with increasing h. 


D. Interdot coupling and quantum phase transition 

Finally we consider the situation where the dots are 
also directly tunnel-coupled, v 0. This generates a 
direct inter-dot exchange J^ir ^ /U which is antiferro¬ 
magnetic. Together with the lead-mediated RKKY cou¬ 
pling, the dots feel a mutual exchange Jtot = TRKKY+-^dir 
which can be tuned from ferromagnetic to antiferromag¬ 
netic by increasing the interdot tunnel-coupling v. The 
result of this is that a high-spin PQD state is realized 
for V < Vc (corresponding Jtot < 0), while a low-spin 
PQD state forms for v > Vc (Jtot > 0). These distinct 
phases are separated by a (level-crossing) quantum phase 
transition, arising at Vc for which Jtot = 0. 






























11 



FIG. 9. Comparison of impurity-DMFT results (red dashed 
lines) and exact NRG (solid black lines) for the PQD model 
with N — 2 (left panels) and 3 (right panels) with hnite in¬ 
terdot coupling u / 0 across the quantum phase transition. 
Upper panels: high-spin phase with v — 0.00511 < Vc. Lower 
panels: low-spin phase with v = O.OID > Vc- Shown for hxed 
U/D = 0.5, e = -U/2 and V/D = 0.06 (but h = 0). 


Interestingly, the Kondo physics for even and odd num¬ 
bers of dots, N, is therefore very different, as discussed 
below. In Fig. 9 we consider specifically N = 2 (left pan¬ 
els) and N = 3 (right panels), in both the high-spin phase 
V < Vc (upper panels) and the low-spin phase v > Vc 
(lower panels). Exact NRG results for the T = 0 spec¬ 
tral function are presented as the solid black lines. 

In the high-spin phase v < Vc (encompassing also the 
case = 0 discussed in Sec. V A) the effective inter-dot 
ferromagnetic coupling Jtot < 0 leads to a combined PQD 
dot spm-N/2 which is Kondo underscreened to spin-(A/' — 
l)/2 below Tk- This holds for even and odd N. The 
inclusion of finite v in this phase does not change the 
spectral lineshape, which is universal for all |ci;| <C |Ttot| 
(only the value of Jtot is changed). In particular, the 
low-energy spectral pinning condition Eq. 19 still holds. 

Including finite v in the impurity-DMET solution does 
not change the results qualitatively, although we note 
that the Eermi level value of the spectrum spuriously di¬ 
minishes with increasing v — see upper panels of Eig. 9. 

By contrast, in the low-spin phase v > Vc where 
Ttot > 0 is antiferromagnetic, even and odd N behave 
differently. Eor even N, the PQD forms a combined 
spin-singlet, and decouples from the leads on the scale 
of |Ttot|- As such, there is a total (on-the-spot) collapse 
of the Kondo resonance on entering the low-spin phase, 
Im FiocGe^^(<^ = 0) = 0. This is seen from the exact 
NRG results for = 2 in the bottom-left panel of Eig. 9. 
Eor odd N, the PQD forms an effective doublet state on 
the scale of iJtotl? which is exactly screened by a stan¬ 
dard Kondo effect below Tk • This leads to a pronounced 
spectral resonance, as seen for A" = 3 in the bottom right- 
panel of Eig. 9. In this case —Im FiocGg^P(O) < 1 because 


the model is not particle-hole symmetric for finite v. An 
effective ESR^^ holds, —Im FiocGe^^(O) = sin^(|nf^p), 
where here the ‘excess charge’of the PQD even orbital 
^fmp 7 ^ f depends on v. 

Impurity-DMET fails qualitatively to capture this 
quantum phase transition, as demonstrated by compari¬ 
son to the red dashed lines. On increasing there is sim¬ 
ply a smooth crossover, with Im riocGg'^p(cj — 0 ) ^ 0 
8iS V ^ oo. In particular, the effective single-impurity 
Kondo physics for odd N and is totally ab¬ 

sent within impurity-DMET. In this case, non-local self¬ 
energy contributions control the physics, and impurity- 
DMET inevitably fails. 


VI. CONCLUSION: APPLICABILITY OF 
INHOMOGENEOUS DMFT 

Although the self-energy becomes strictly local only in 
the limit of infinite dimensions,^ DMFT employing local 
self-energies is widely used in the approximate treatment 
real systems with finite coordination.^’^ Furthermore, in¬ 
homogeneous (or real-space) DMFT has been applied in 
diverse contexts, including e.g. layered materials^ and 
molecular electronics. 

DMFT is certainly a very powerful tool, providing deep 
insights into the physics of complex systems.^ However, 
for highly inhomogeneous systems where both strong in¬ 
teractions and details of real-space geometry are impor¬ 
tant, it is often not obvious whether a local self-energy 
approximation could yield physically sensible results. Of 
course, the problems to which inhomogeneous DMFT is 
applied are typically those for which exact results are 
not available, and so the quality of the approximation 
can remain unclear. 

In this paper we have applied “impurity-DMFT” to 
the 2 IAM (featuring two magnetic impurities separated 
in real-space on the surface of a 3D cubic lattice), and 
PQD models (in which quantum dots are coupled in par¬ 
allel between source and drain metallic leads). The local 
self-energy, combined with self-consistency, constitutes 
a rather sophisticated approximation, and nontrivial re¬ 
sults are obtained. 

Although the 21 AM and PQD models are among the 
simplest examples for which DMFT can be used, their 
exact solution with NRG is still very challenging.^^ By 
comparing exact NRG and impurity-DMFT results for 
dynamical quantities in these models, we were able to 
test the validity of the local self-energy approximation in 
highly inhomogeneous contexts. 

We find that in situations where Kondo screening is 
the dominant effect, impurity-DMFT can work surpris¬ 
ingly well, especially at low energies. For example, in 
the 21 AM, impurity-DMFT correctly predicts the onset 
of the dilute limit,and recovers the Fermi level value of 
the impurity spectral function. In PQD models, subtle 
underscreened Kondo physics is apparently captured in 
impurity-DMFT through the self-consistent generation 
of a critical pseudogapped effective medium. Impurity- 
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DMFT is also correctly able to describe the effects of an 
applied magnetic field in PQD systems. 

However, we also uncover cases where the local self¬ 
energy approximation fails catastrophically. These are 
situations where inter-impurity correlations are responsi¬ 
ble for suppressing Kondo physics. For example, effective 
antiferromagnetic RKKY interactions in the 2IAM can 
lead to inter-impurity singlet states which simply cannot 
be described within impurity-DMFT. In PQD models, 
the quantum phase transition between high- and low- 
spin states on tuning inter-dot coupling is similarly not 
reproducible in impurity-DMFT. These results naturally 
have implications for the wider applicability of DMFT to 
inhomogeneous systems. 

Finally, we comment that exactly-solvable quantum 
impurity problems provide an ideal context to develop 


and test novel extensions of DMFT, going beyond the 
local self-energy approximation. Reliable solutions will 
then be within reach for a wider range of problem. 
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